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Abstract 

A six-dimensional Navier-Stokes equation derived by one of the authors(ST) is solved 
for a turbulent mixing layer to demonstrate that it has a solitary wave solution. Tur- 
bulence intensities and Reynolds' stress are calculated using this solution, showing satis- 
factory agreement with experiments although no emperical constants are involved in the 
theory. 

Key words: turbulence, solitary wave, mixing shear layers 

With their recent advent parallel/superparallel computers have provided a power- 
ful tool for computational fluid dynamicists. With regards to direct numerical simula- 
tion(DNS) of turbulence, however, the gap between the Reynolds number 10 4 that can 
be reached using those computers is still far short of the target ones 10 7 ~ 10 8 to be 
practical with airplane designers. This computational difficulty reflects the fact that with 
increase in the Reynolds number R, grid size required for DNS becomes finer as R~ 3 ^. 
This situation is caused by the fractal structure of turbulence ]TJ0| that the Kolmogorov 
scale(~ i? -3//4 ) is the lower fractal limit above which turbulent fluid quantities are not 
considered as smooth, but are self-similarly irregular with fractal dimension of about 2.36. 
This is where differentials are not to be replaced with finite differences, or desperately 
large number of Fourier modes need to be considered in employing the spectral method. 
To remedy this 'small eddy difficulty' of high Reynolds number turbulence, it is advisable 
to work with averaged fluid variables rather than fractal ones directly. Of two candi- 
date methods to meet this purpose, namely, the large-eddy simulation and the statistical 
theory of turbulence, the latter which is of classical origin || is reexamined on the basis 
of first principles of non-equilibrium statistical mechanics^] ~ M: It has led us to an 
averaged Navier-Stokes equation in a six-dimensional(6D) space comprising 3D physical 
plus 3D eddy spaces ||. 
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This formalism allows us to deal with smoothed quantities at the expense of increased 
number in independent variables. Furthermore it has a computationally favorable feature 
that the small eddies are mapped onto a small region around the origin of the eddy 
coordinates. Thus the small eddy difficulty is eliminated by employing fine mesh only 
in the vicinity of the origin of the eddy space, not necessarily the whole computational 
domain as in the DNS. 

The objective of this paper is to check if the equation thus derived has a solution that 
is physically sound by solving a problem of practical interest. A turbulent mixing shear 
layer serves as such an example. 

Since the discovery of the coherent structure of turbulence [ID ], this problem has at- 



tracted attentions of computational fluid dynamicists as well as experimenters. DNS, 
however, has not been able to reach Reynolds numbers high enough to establish spatial 
self-similarity. Efforts are then directed to solving a similar flows of tangentially opposing 
parallel streams where turbulent mixing layer grows timewise, not streamwise. Reliable 
DNS result is reported up to R = 2 x 10 4 (JTTJ . This flow configuration is much easier to 
compute(average flow does not change streamwise). This parallel mixing layer, however, 
is fictitious in nature(no experiments possible), and the similarity to the spatially growing 
mixing layer is only qualitative(no coordinate transformation of equations and boundary 
conditions from one to the other is possible). Therefore for quantitative comparison with 
existing data to be possible, we need to solve spatially growing case directly. 

Experimentally a turbulent mixing layer is generated between two parallel streams 
with different velocities U±oo, starting to merge at X\ = 22 = 0. It is shown [[L2[][]13] that 
2-D incompressible mixing layers tend to be self-similar with increase in thickness 5 of the 
mixing region that grows linearly with streamwise distance; 5 = ax\. This experimental 
evidence indicates that the governing equations [Eqs.(34) and (35) of ref.9] 

d j Q j = (1) 
(j (j i f)a ■ 

- c ^ + u ^ - vd t)Qj + + ir-Q* + = ( 2 ) 

Ot Ooi p OXi 

where Q replaces p~ x q(p ; the density) of ref.9, and 

are controlled by reduced number of self-similar variables 

V = x 2 /5, s = S/5 (4) 

where subscripts 2 and 3 refer to normal and spanwise directions viewed from the splitter 
plate separating the two parallel streams. This observation of self-similarity together with 



the group-theoretical consideration [14fl tell us that terms with viscosity be negligible, 
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corresponding to the physical situation that the viscous stress is overwhelmed by the 
Reynolds stress everywhere in the flow. 

Even under this simplifying condition, the equations have four independent variables 
to which existing computational methods are not applicable directly. To make them 
tractable it is proposed to suppress one of the eddy variables; d/ds2 = 0. This choice 
saves the minimum requisite that longitudinal vortices depend on S3 and transverse ones 
on Si more strongly than others. Under these conditions Eqs.(l) and (2) reduce to the 
following set of 'inviscid' equations: 

a . dq 2 dq 3 

NL qi + d±q4 — arju'qi + u'q 2 = (6) 
d 

NL q 2 + —q A - ar]v'q 1 + v'q 2 = (7) 
or] 

NLq3 + ^ = (8) 
Os 3 



with 



d d d 

di = (1 - as^- oc(rj— + s 3 —) 

os 1 or) 0S3 

NL = — - c— + {u + qi )d! + {v + q 2 )— + q 3 — 
at asi or) os 3 



(9) 



where (u(r)),v(r)),0) is the mean velocity vector made nondimensional using the velocity 
difference U = f/oo — f^-oo ,and q a and c are nondimensional fluctuations and phase velocity 
defined similarly; 



qj = Qj/U q 4 = Qi/pU 2 
c = C/U . 



(10) 



Since this paper is focused on a quick check of the adequacy of the proposed approach, 
in particular, the existence of solitary wave solution, we do not solve, but instead adopt 
measured values for u, while v is computed using the continuity equation; 

n du 

v = r)—ar). (11) 
Jo ar) 

The phase velocity c, namely, the propagation velocity of turbulent vortices in the 
(x, s) space is yet to be specified. Whether it is to be given on first principles is not well 
understood at this stage. Therefore we will explore this point by invoking two different 
assumptions representing two extreme cases of large and small eddies, and will see how 
they compare to each other. 
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The one is to put 

C/oo + U- c 



- (12) 
whose right-hand-side represents the phase velocity of vortices as observed by visualisation 
experiments, in other words, the propagation velocity of large eddies The other is to 
put 

c = u (13) 

which represents Taylor's hypothesis that claims the vortices to be carried with local mean 
flow. It is more likely that small eddies obey (13) rather than (12). 

Initial-value problems are solved employing the pseudo-compressibility method[l6|] in 
the 3-D space (r),si,s 3 ). The method assumes a fictitious term (5~ x df^jdt {(3 3> 1) 
added on the continuity equation (5) to allow the whole set of equations (5) through (8) 
for time evolution procedure. 

A unique solitary wave solution is shown to build up eventually independent of what- 
ever initial distributions to start with. The solution having the property 

51,92,54 ; even in s 3 , 
q 3 ; odd in s 3 

is sought to meet with the two-dimensional structure of the mean flow; (u'iU 3 ) = (u' 2 u 3 ) = 
0. It is also confirmed that the time-dependent method leads asymptotically to a steady- 
state solution, in other words, that there are no signs for the solution to have a limit-cycle 
or strange-attractors. 

The mean flow condition employed for the computation is 

U-oo/U^ = 0.6 (14) 



the same as the experimental condition of refs.[[L2(] and [T3|]. The mean velocity profile for 
u is taken from ref.|l2|, and that for v is calculated from Eq.(ll). 



In Figs.l are shown the birds-eye- view representation of the solitary wave solution 
for 5i(?7, Si,s 3 ) projected onto the plane of r) — 0. The computation is based on two 
different assumptions(12)(Figs.la) and(13)(Figs.lb) for the phase velocity. They look 
quite similar at least qualitatively. In both streamwise(si) and spanwise(s 3 ) directions of 
eddy coordinates, the main peak is accompanied by valleys, which is the sign of existence 
of adjacent counter- rotating vortices. 

Once the solitary-wave solution has been obtained the fluctuation correlations and the 
turbulence intensities including the Reynolds stress are calculated through the following 
integrals 

(u'j^ueix)'} = Rji(x,x) 

1 f f°° Xi — Xi x 3 — x 3 . 



(2n) 



f f°° /- X\—X\ x 3 — x 3 
J J ds 1 ds 3 q j {r],si,s 3 )q e {r],s 1 ^ , s 3 H ) (15) 
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which is the 2-D and self-similar version of formula(39) of ref.9. Figs. 2 show the turbulence 
instensities obtained by putting x = x in (15): Those figures correspond to a) streamwise 
(j = £ = 1), b) transverse (j = £ = 2), c) spanwise (j — £ — 3) intensities of velocity 
fluctuations and d) Reynolds stress (j = 1, £ = 2). Two curves with solid lines correspond 
to respective assumptions for the phase velocities as indicated in the figures. It is observed 
that both theoretical curves fit better with Bell-Mehta's data rather than Wygnanski's. 
For each case they are reasonably close to each other considering that the two assumptions 
stand for the opposite extreme of the eddy size. This computational evidence implies 
robust tendency of the governing equations towards the solitary wave formation. As is 
easily confirmed, the governing equations do not allow for solutions that are symmetric 
with respect to rj. It may well be expected that the seemingly symmetrical distribution 
as observed in experiments makes its appearance only after the integration over s-space 
has been effected (formula (15)). The sensible asymmetry recognized in computed results 
may be attributable to insufficient mesh numbers as limited by the memory size of the 
computer. 

Fig. 3 shows the vertical turbulent transport of the Reynolds stress #122(77) =< u' x u' 2 U2 > 
for c = Uq and compared with experiment [Q|. The agreement between theory and exper- 
iment for such a small and subtle quantity is unexpectedly close to each other. 

Fluctuation-correlations of the vertical velocity component R22 = (u' 2 u' 2 ) using (15) in 
streamwise (a;r/ = 0, 23 = £3) and spanwise (b;^ = 0,xi = xi) directions, respectively, are 
computed and plotted in Figs. 4. Each of these figures corresponding to the two postulates, 
c = Mo and c = u, falls very close to each other, so only the case c = u is shown here. 
Both correlation curves have the points of zero-crossings, indicative of adjacent vortices 
of opposite rotations, also of the clearcut existence of bulk of co-rotating cores in both 
directions. Length £ 3 of zero-crossing in the spanwise correlation curve is available from 



Bell-Mehta's experiment [O] , falling quite close to the predicted one. 
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Figure Captions 



Figure 1: Solitary wave solution of qi(r), s±, s 3 ) on two different assumptions for phase 
velocity c [(a) and (b)] viewed on the plane of rj = 0. 



Figure 2: Distributions of turbulent intensities in transverse (77 -) direction, compared 
with data by Oster-Wygnanski[12] and Bell-Mehta[13]; a)streamwise (Ru), b)transverse 
(R22), c)spanwise (-R33) and d)Reynolds' shearing stress (i?i2)- 



Figure 3: Transverse transport of Reynolds' stress (-R122) and comparison with data 
by Bell-Mehta[13]. 



Figure 4: Stream- and span-wise variation of correlation of the vertical velocity com- 
ponent R22 at i] = 0. 
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